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: 

Transport properties of narrow two-dimensional conducting wires in which the electron scatter- 
ing is caused by side edges' roughness have been studied. The method for calculating dynamic 
characteristics of such conductors is proposed which is based on the two-scale representation of the 
mode wave functions at weak scattering. With this method, fundamentally different by-height and 
by-slope scattering mechanisms associated with edge roughness are discriminated. The results for 
single-mode systems, previously obtained by conventional methods, are proven to correspond to the 
' former mechanism only. Yet the commonly ignored by-slope scattering is more likely dominant. 

, The electron extinction lengths relevant to this scattering differ substantially in functional structure 

^ . from those pertinent to the by-height scattering. The transmittance of ultra-quantum wires is calcu- 

^ ' lated over all range of scattering parameters, from ballistic to localized transport of quasi-particles. 

The obtained dependence of scattering lengths on the disorder parameters is valid qualitatively for 
arbitrary inter-correlation of the boundaries' defects. 
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PACS numbers: 72.10.-d; 72.15.Rn; 73.20.Fz; 73.23.-b 



INTRODUCTION 



o , 

^ O ^ , The influence of inhomogeneities and defects of different types on transmission properties of waveguides for both 
quantum and classical waves is a subject of long-term intensive research. Most of theoretical studies have been 
restricted to flat surfaces. In practice, however, side boundaries of mesoscopic conductors are either inherently rough 
^ (due to growth, fracture, etc.), or artificially patterned (e.g. due to lithographic preparation). One way or the other, 
surface-corrugated conducting systems find widespread applications in technology and material science (quantum 
conductors, optical fibers, etc.), and, hence, it is often necessary to know how the shape of the confining surface 
, affects the charge or classical wave transport in different systems. 

' Of peculiar interest for contemporary microelectronics, from both applied and theoretical viewpoints, are two- 
. dimensional (2D) conductors of mesoscopic size. Nowadays there exists widely established opinion (see, e.g., Refs. 
H, ||, ^ 1^, ^, 0, H 1^, |l^, |ri|, ^ and references therein) that dynamic properties of pure-in-bulk quantum wires, in 
particular two-dimensional, are largely determined by scattering the electrons from randomly rough side boundaries 
of the conductor. This scattering mechanism is proven to be responsible for both relaxation processes in multi-mode 
conductors ||^, |^, ^ and for non-dissipative (Anderson) localization of conduction electrons in narrow single-mode 
wires §0. 

^ , In studying the electron transport in 2D surface-corrugated systems two main problems are especially highlighted. 
O ■ One of them is relevant to adequate description of the electron scattering from statistically rough surfaces. The other, 
dynamic, problem is pertinent to proper consideration of the interference of multiply scattered quantum waves, which 
is essential for describing within the framework of a perturbation theory the effects resulting from non-dissipative 
localization of the electron states. 
}J] ' To resolve the first problem and make an adequate comparison with experimental data one needs a theory which 
5^ ' relates transport properties of the conduction electrons to the shape of bounding surfaces of the conductor. Since 
in random-inhomogeneous 2D wires fluctuations of both of the side boundaries can be considered as either mutually 
independent or correlative, subject to the preparation technology, it is important to trace the relation between the 
kinetic quantities and the statistics of boundary irregularities, as well as mutual correlation of the opposite boundaries. 



In the previous paper 1 10 we have studied the case of a single- mode 2D conductor with statistically identical rough 
side boundaries. The particular model considered in Ref. of the wire with completely correlated boundaries 
(CCB) is equivalent to the deterministic waveguide system of constant width whose inhomogeneities consist in solely 
waveguide bends. It was shown that the electron dynamics in such conducting strips is governed by quite a different 
Hamiltonian than that pertinent to the seemingly more general model of 2D conductor with one boundary which is 
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rough and the other being ideally flat ||^ . A substantial distinction of the Hamiltonians in Ref. and Ref. [|o) have 
resulted in qualitative functional distinction of the obtained scattering lengths. 

The results of Refs. and have provided the grounds to propose that in those works the quantum wave 
scattering can be assumed to be associated with fundamentally different physical factors, namely with deviation 
of the boundaries from their "ideal" shape in Ref. and fluctuation of their slopes in Ref. [|lOj. Provided the 
supposition is true, when studying particle or classical wave transport in surface-corrugated waveguide systems one 
has to distinguish between two different non-interfering scattering mechanisms which we call by-height (BH) and by- 
slope (BS) scattering. However, it should be noted that quantum wave scattering was analyzed in Refs. and |]lOj for 
different waveguide geometries and with the use of substantially different methods. Therefore, the conjecture stated 
in Ref. about the relative significance of these scattering mechanisms needs to be additionally substantiated. 

To avoid possible misunderstanding and support our idea of different scattering mechanisms pertaining to the 
imperfect boundaries of guiding systems, we examine the waveguide (conductor) geometry admitting of both BH and 
BS scattering simultaneously. We consider a 2D conducting strip with statistically symmetrical rough boundaries (the 
abbreviation SSB will be used for such a strip, in contrast to the CCB strip considered in Ref. 0), which is physically 
equivalent to the waveguide with a straight central line (guiding axis) and a randomly fluctuating width. It will be 
shown below that the suggested model of the waveguide corresponds qualitatively to a 2D wire with opposite side 
boundaries whose inter-correlation can be thought of as arbitrary. For the SSB model, as well as for the arbitrary 2D 
waveguide, the presence is typical of both BH and BS scattering mechanisms. It is noteworthy that these mechanisms 
compete with one another, depending on the roughness statistical parameters, even in the simplest case of boundary 
asperities being small in height and rather smooth. 

Technically, specifying the scattering mechanisms is most straightforward by reducing the problem of the electron 
scattering from complicated boundaries of the conductor to the appropriate "bulk" problem specified by a complex 
Hamiltonian but simple boundary conditions. Analysis of the problem thus formulated can often be found much 
easier than of the problem with complicated boundary conditions. In some cases it can even be performed non- 
perturbatively. In this work, using such an approach, we managed to reasonably discriminate between BH and BS 
scattering mechanisms and analyze their competition in the electron-surface scattering. Besides, solving the "surface" 
problem in "bulk" formulation enabled us to carefully trace such a fine spectral effect as the Anderson localization of 
current carriers. 

Note that the method for solving the problems of the by-surface scattering through reduction to the Hamiltonian 
form is not quite original. It was employed, in particular, in Refs. and j^] where the coordinate transformation was 
used, which smooths out the rough surface to the flat one. In this paper, analogous reduction of a "surface" problem 
to the Hamiltonian formulation is made by merely going over to the local mode representation. This approach ensures 
the optimal choice of trial quantum states which serve as a basis for the perturbation theory. In our view, it is essential 
that the mode states a priori contain the information on lateral confinement of the system under consideration and, 
therefore, are more adjusted to perturbative treatment of transport problems in waveguide-like systems than the 
widely used isotropic plane- wave basis. 



STATEMENT OF THE PROBLEM: CHOOSING STATISTICAL MODEL 

We consider a conducting strip of average width D with the non-uniform stretch of length L (shaded region in 
Fig. 1) occupying in {x,y) plane the area restricted by the inequalities 

- L/2 < X < L/2 , 

(1) 

-D/2 + ^,{x) < y <D/2 + U^) . 

Functions ^1,2 (a^) describing the side boundaries' roughness are assumed to be random processes with zero mean 
values. It turns out to be more convenient for subsequent analysis to introduce, instead of ^1,2(2^), two different 
random functions, viz. 

w{x) = D + Aax) and U^) = [^,{x) + ^2{x)]/2 , (2) 

where A^(a::) — $2(2;) — ^i(x). Function 'w{x) is nothing but the fluctuating width of the strip whereas S,c{x) describes 
the y coordinate of the randomly fluctuating central (symmetry) line of the 2D waveguide. In terms of the functions 
(^, the disordered region (|l|) is represented by inequalities 

- L/2 < X < L/2 , 

(3) 

-w{x)/2 + S,c{x) < y <w{x)/2 + i^{x) . 



FIG. 1: 2D electron waveguide with rough side boundaries pertaining to the SSB class. 



Reformulation of the problem in terms of new characteristic functions is advantageous from the mathematical point 
of view. Besides, it gives a clear idea of the part played by different physical factors in quantum wave scattering in 
surface-corrugated conductors. 

Bearing in mind the statistical nature of the problem under consideration it is necessary to specify correlation 
properties of new random variables in conformity with the conductor geometry chosen. The mean values of 
the functions ^/{x) and S,c{x) are naturally {w{x)) = D and {^dx)) = 0, the angular brackets stand for statistical 
averaging over realizations of ^1,2(2^)- Unlike the plain averages, the binary correlators of functions (|^) are not so 
uniquely predetermined but depend significantly on the intercorrelation between the side boundaries of the conducting 
strip. Here we consider the correlation model in whose frames both of the boundaries are regarded statistically identical 
in that the functions ^1.2(2^) obey the equalities 

{^^{x))=0 and {^,ix)^,ix')) =a^Wix-x') , i ^ 1,2. (4) 

Here a is the rms height of the asperities, which is thought of to be identical for both of the strip edges, W(a;) is the 
correlation coefficient specified by the unit maximal value and the correlation radius Tc- 

As for the intercorrelation of the opposite-boundary roughness, two marginal options seem to be distinct among 
others for their particular symmetries. One of them, abbreviated in Ref. as CCB model, in terms of statistical 
formulation implies the fulfilment of correlation equality {^i{x)^k{x')) = a^yV{x — x') at z ^ fc; it is physically 
equivalent to a 2D waveguide of constant width, whose inhomogeneities are in the form of random bends. 

In this work, we make use of another particular model in which the displacements of the side boundaries comply 
with the first of relations (Q) whereas for the opposite boundaries the correlation equality holds 

{x)(.k {x')) = -n^W{x - x') , i^k. (5) 

It is easy to make sure that equations (^) and ^ are equivalent to the following correlations, 

{Ux)w{x')) = , (6a) 

{Ux)Ux'))^0. (6b) 

The equality ( |6a| ) implies that the functions £,c{x) and w{x) can be thought of, within the correlation approximation, 
as statistically independent random processes. The second equality, (|6b[), within the same approximation is consistent 
with deterministic relation ^i(x) = —^2(2;), that holds true in the case of the electron waveguide with side boundaries 
fluctuating symmetrically about the straight central line, as is shown in Fig. 1 (the SSB model). 

From the linear response theorypTt, the dimensionless conductance g{L) (in units of /ttK) at T = is represented 
by the expression 



,T^ ^ ft ^ ^ ,dG{ry)dG*{r,r') 



Here G{r, r') is the retarded one-electron Green function; integration with respect to r = [x, y) runs over the area (|^) 
occupied by the irregular part of the wire. Within the isotropic Fermi liquid model, the function G(r, r') is governed 



by the equation 
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kl+iO) G{r,r') = S{r-r') 



(8) 



where A is the two-dimensional Laplace operator, kp is the electron Fermi wave number. As to the boundary 
conditions to equation (|^), in the direction of current (x) the conductor will be regarded as open whereas in the 
transverse direction (y) the zero (Dirichlet) conditions will be used. Note that formula ( [7|) i s valid asymptotically 
under conditions of weak scattering (WS) [ [T8| , [Tot , which we assume in this paper (see Eq. (31) below). 

Inhomogeneity of side boundaries of a waveguide system can be taken into account in a number of ways. The most 
commonly used method reduces to linearization of the exact boundary conditions in small fluctuations of bounding 
surfaces p^, pO| . However, such an approach can result, in some cases, in significant underestimation of the inter-mode 
scattering which often exerts primary control over the dynamics of waves within bounded regions. It will be shown 
below that the entanglement of waveguide modes is governed by fluctuation of the confining surfaces in slope rather 
than in their displacement from the specified ideal shape. Meanwhile, if one linearizes the boundary conditions, 
it becomes precisely the height of the surface roughness that serves as a determining factor for wave scattering. 
Therefore, entanglement of the modes can be, either to a large extent or even entirely, lost when expanding the 
boundary conditions in small heights of the surface roughness. 

The inter-mode scattering can be taken into account using the methods of Refs. jlj, |l^ or with the method of 
"smoothing" coordinate transformation, used in Refs. [0, §0. In the latter papers, the dynamic problem pertinent 
to the system with corrugated boundaries is reduced without any approximation to the analogous problem related to 
the system with ideal boundaries, though governed by a more complex Hamiltonian. This method is advantageous 
in that it enables one to analyze the dynamics of quasi-particles without resorting to the concept of "adiabaticity" of 
the confining potential |^ |2| . 

Technically, the procedure of "smoothing" rough boundaries of the conductor is most convenient to perform without 
applying explicitly the coordinate transformation, which, in addition, is quite non-linear in general case. The same 
result can be obtained by going to the local (in lengthwise coordinate x) mode representation in Eqs. (^ and (||), 
i.e. by performing Fourier transformation in the coordinate y using the complete set of "transverse" eigenfunctions 
of the Laplace operator which are consistent with boundary conditions prescribed at the true walls of the wire, 
y = ±w{x) + $,c{x). In our case (Dirichlet conditions) we choose, for definiteness, the functions 
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By substituting G(r, r') in the form of a double Fourier series, 



G(r,r') 



n.n' — 1 



Sn{y\x)Gnn'ix,x')Sn'{y'\x') , 



(9) 



(10) 



into equation (g|) we arrive at the following set of equations for the coefficients G„„' of that series which we call 
hereafter the mode Green functions. 
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w{x) 
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(11) 



The numerical coefficients in equation {ulh have the form 
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l^(,w{x) are the differential operators of the type 



d d 

U^{x) =w\x)— + —w'{x) 



(13a) 
(13b) 



The conductance expression (|7|), on substituting there Green function ( |lO| ) and integrating over coordinate y, is 
reduced to the form 



J-L/Z J-L/2 n,n' = l 



dx 



E 

m = l 



dx' 



E 

m' = l 



Gl^,{x,x')<^>„^'n'{x') 



(14) 



Here A'^c = [kD/ir] is the number of "open conducting channels", that is extended waveguide modes. The coefficient 
matrix {^nm} in (0) is composed of the elements 



^nrn{x) = 2B„ 



u/(aO 

I / \ V-'- 

w(a;) 



i / \ 1 

W[X) 



(15) 



whence it is evident, subject to Eq. ([l^), that there is no diagonal terms, <i>„„(x) = 0. 

We emphasize that the set of mode equations (|ll|)-(|l5|) is applicable for any correlation model of rough boundaries. 
Equation (^l|) for the mode Green functions contains the effective electron-surface scattering potential which consists 
of two (substantially different in their physical meaning) types of terms. The first term, [7m/w(a;)]^, is determined by 
random local displacement of the waveguide boundaries and, therefore, is responsible for the by-height (BH) scattering. 
All the other terms contain the gradients either of the wire width, w'{x), or of the y-coordinate of the symmetry line, 
Cci^)- Consequently, they describe the by-slope (BS) electron scattering. Note that the structure of the BS part of 
the scattering operator is not highly sensitive to the difference between two specific models of the waveguide, when 
either w'{x) = or ^ ^i^) = 0. At the same time, depending on these models, the role of the BH term in the electron 
scattering changes drastically. Indeed, if w'{x) = (i.e., w{x) = const) the electron scattering arises only due to 
fluctuations of the roughness slopes that are described by the BS terms with ^'^(a;). Otherwise, when Cd^) = 
both the BS and BH scattering potentials are contributing. The case of w'{x) = corresponds to the CCB model of 
the waveguide |jl^ (constant width, bend-type inhomogeneities) whereas the case Cd^) = is physically equivalent 
to the SSB model with the correlation properties (|4|)-(||) ("straight" waveguide with fluctuating width). So, we can 
conclude that the SSB model is qualitatively similar to the general one with arbitrary correlation properties of the 
lateral boundaries. On the other hand, employing this model allows one to avoid excessively cumbersome calculations 
without sacrificing the quality of the results. 



REDUCTION TO ONE-DIMENSIONAL DYNAMIC PROBLEM 



In studying the electron transport in bounded systems, one of the basic characteristics is the number of conducting 
channels or, what is the same, the extended waveguide modes. However, the number of eigen-modes cannot in all 
cases serve as a good quantum parameter if one addresses the waveguide with a variable cross-section. Nevertheless, 
formulation of the investigated problem in mode representation is advantageous, since in this approach the mode 
structure of the conductor can always be fixed, = const, while boundary displacement can be considered as a 
source of perturbation of the mode Hamiltonian in Eq. (^l]) . 

In this work, bearing in mind the perspective of Anderson localization of the current carriers owing to such "edge 
perturbations" , we consider a single-mode 2D strip with the average width D confined within the interval 

TT/kp < D < 27r/fcF ■ (16) 

For the reasons being discussed below (see also the Appendix), fluctuations of the conductor width will be regarded 
small as compared to its average value. 
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Under restrictions (|T^), the only element of the mode Green function matrix {G„n'} whose contribution to the 
conductance at weak scattering (see Eq. (^T|)) is not parametrically small is the intra-mode propagator Gii{x, x'). 
Before deriving the closed equation for this function, it is worth noting that under certain conditions the operator 
potentials U^^j^{x) can be left alone in square brackets standing in front of the inter-mode propagator Gmn' in equation 
(pl|). Indeed, the relative part of the terms quadratic in roughness slope and the linear terms in those brackets is 
estimated by the ratio 



wix)\k^'{x)+C'ix)\ D{l + krc) ■ 
It is evident that if the rms height of the boundary displacement is small enough, 

cj/D < 1 , 



(17) 



(18) 



one can neglect the terms that are quadratic in ^' in comparison with their adjacent linear counterparts. Besides, for 
the sake of simplicity of the subsequent calculation we will consider the boundary asperities to be also smooth. 



cr/r-c < 1 , 



(19) 



yet without requiring their adiabaticity in the sense of Refs. |2^ . 

The limitations and (|l|) are common in solving the problems of classical wave scattering using the perturbation 
theory |]l3|. The restrictions are motivated by the necessity of eliminating the well-known "shadowing" effect. In view 
of likeness of the mathematical technique, this issue arises in the considered quantum problem as well. However, 
under the conditions ( [l8| ) and (|l^) there does not exist the wave shadowing in a single-mode waveguide. 

With all the above reasonings taken into account, we obtain the equation for Green function Gii{x,x') setting 
71 = n' = 1 in Eq. (0), 



dx'^ 



+ kl+tO-Vhix)-Vsix) 



Gll{x, x') - ^ Ulm{x)Gml{x, x) = 5{x - x) . 



(20) 



m=2 



Here we have introduced the notation kf for the "unperturbed" lengthwise energy of the extended mode n = 1, 



kl 



hp 




(21) 



The quantities Vh{x) and Vs{x) in square brackets of Eq. (|20|) stand for the intra-mode potentials and have the form 



w'^{x) 



Vh{x) 



Vs{x) = (1 + y 




w'{x) 
2w{x) 



(22a) 
(22b) 



Assuming these potentials to be taken as a perturbation, we design them purposefully to have {Vh,s{x)) = 0. As to 
the inter-mode operator potential in Eq. (EQ), 



Ulra{x) = -Bir 



w{x) 



Uw{x) 



(23) 



it possesses this property by definition. 

Although the mode Green functions in Eq. (|^ depend on a single space variable, the problem certainly cannot 
be thought of as one-dimensional for coupling the function Gii{x, x') to all of the inter-mode propagators with mode 
indices m 1. Nevertheless, one can obtain a closed equation for this function, which at weak scattering admits 
of the perturbation analysis adequate to the diagrammatic method of Refs. ||2^ and ||l^. In Ref. [^, the evidence 
was given that in the case of arbitrary 2D imperfect waveguide system governed by a set of dynamic equations 
of the same functional structure as that of the system (pT|), all of the inter-mode propagators Gmnix,x') can be 
expressed, by means of a linear operator, in terms of only one intra-mode Green function Gnn{x, x'). Although 
for arbitrarily disordered systems this relation is of little practical value for the complex dependence on scattering 
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potentials, under weak scattering conditions the reciprocal operator expressions prove to be rather uncomplicated and 
readily analyzable. Referring an interested reader to Ref. psf for the exact procedure, we give here a simple recipe 
for obtaining the approximate relation between the inter- and intra-mode propagators, which is valid in the case of 
weak electron scattering. 

This recipe, already used in Ref. , consists in solving the set of equations (|l^) iteratively with respect to the 
inter- mode propagators entering equation (|2^). By letting n' = 1 and re-designating the mode variables one can 
reduce Eq. (nTl) to a set of inhomogeneous equations with respect to the functions G„ii{x,x') with m > 1, 



g„i^G„ii{x,x') - ^ U„-,k{x)Gki{x,x') ^ Ujni{x)Gii{x,x') 



(24) 



k = 2 



It can be easily seen that all interesting propagators Gmi can be expressed linearly through the single intra-mode 
propagator Gu. In equation (p^, is the differential operator from curly brackets of Eq. (pi]), where the term 
with ^c^(a;) should be omitted, in accordance with the SSB model chosen, and replacement of mode indices should be 
made n ^ m. The inter- mode potentials Umk{x) in Eq. ( p^ ) have the form similar to the operator (p3|). 



U,nk{x) = -B„ik—r^iiw{x) . 

W[X) 



(25) 



It is essential that in the case of weak scattering all Green functions Qrn{x,x'), with the aid of which the equation 
(^) is to be solved, belong to the class of so-called evanescent functions. They are sufficient to be taken in the 
unperturbed form (i.e. to zeroth order in boundary roughness) which are found to be strongly localized. 



Gi^Kx,x') 



■exp(-|fc„i||a;-x'|) 



1/2 



Here |fc„i| = (jTm/D) — kp > 0. Since at weak scattering all of the inter- mode potentials entering Eq. 
be regarded as small (in functional sense), the following approximate relation can be derived iteratively, 

r-L/2 



Gjnl{x,x') 



L/2 



dxig^^'>{x,xi)U,„i{xi)Gii{xi,x') 



By substituting 



dx^ 



into (|20|), we arrive eventually at a closed equation for the intra-mode propagator Gu, 

■L/2 



+ kl + iO-Vhix)-Vsix) 



Gii{x,x') 



dxi K{x, xi)Gii (xi, x') — 6{x — x') 



(26) 
|) can 

(27) 
(28) 



-L/2 



Equation ( Pq ) contains complete information on scattering of single-extended mode electrons by the roughness of 
the conductor boundaries. While local potentials Vfi{x) and Vs{x) are responsible for the intra-mode scattering, the 
integral operator in (ES) accounts also for the inter-mode scattering. The kernel of this operator has the form 



k{x,x') 



w{x) 



u^{x)g^^\x,x')uux') 



w{x) 



UM)S^^\x,x')U^{x') 



(29) 



being, in its turn, a differential operator. Similar to the intra-mode potentials T4,s(x), the operator potential ( p9| ) is 
constructed in such a way as to make (^K{x,x')'^ — 0. This certainly brings about even greater complication of the 
exact form of the longitudinal energy kf in comparison with that given in Eq. (pT|). But the smoothness condition 
( p9| ) makes it possible to omit the terms containing the derivative w'{x) in the "unperturbed" mode energy, thus 
replacing it with the simplified value kf ^ kp — (ii'^ / vJ^ {x)) . 

At last, in formula ( p^ ) for the conductance, in view of the single-mode geometry of the conducting strip, Eq. (p^), 
one should keep only the terms with the diagonal mode propagator Gii{x, x'). With Eqs. ( p^ ) and (p^), the conduc- 
tance expression reduces to a relatively simple form. 



L/2 



dx 



L/2 



, 9Gii(x,x')9G^i(x,x') 

dx 



L/2 



dx 



dx' 



(30) 



which will be the subject of further analysis in conjunction with equation (Eg). 
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TWO-SCALE APPROACH TO THE PROPAGATOR CALCULATION 



Since we focus our attention on the single- mode quantum wires, the unperturbed Green functions (26) of all higher 
modes with m > 2 are real-valued and strongly localized in space. The kernel of the operator potential in equation 
( p8| ) is thus hermitian, which ensures hermicity of the total Hamiltonian of the dynamic system. 

Equation ( p8| ) can be solved asymptotically within WS approximation. This approximation is most convenient 
to be stated in terms of spatial lengths inherent in the problem. We will arrange these lengths into two groups, 
viz. "microscopic" and "macroscopic" lengths. The electron Fermi wavelength, 27rfc^^, and the correlation length of 
boundary asperities, Tc, are assigned to the first group. Among the macroscopic lengths there will be scattering isc of 
the extended mode n = 1 (to be determined below) and the length L of an irregular part of the conductor. In terms of 
these lengths, the criterion for scattering to be classified as weak can be expressed through the following inequalities, 

kp\r,^£sc,L. (31) 

Note that correlation between the lengths pertaining to the same group (either micro- or macroscopic) can be thought 
of as arbitrary. 

The retarded Green function entering the conductance expression ( |30| ) is a solution to the boundary-value problem 
defined by equation (^) and the appropriate boundary conditions. We assume the conducting strip to be open at 
the ends x = ±L/2 in the direction of the current. For an open waveguide system, Sommerfeld's radiation conditions 
well-known in classical wave theory p3| , p5| are appropriate. In the case of ID equation ( p8| ) they can be expressed in 
the Leontovich's form [Ol Wj\ 



d \ 

— =F«fci j Gii{x,x') 



- , (32) 

x=±L/2 

where the source is assumed to be placed inside the waveguide, x' G [—L/2, L/2]. 

For lack of dynamic causality, the solution of Sturm-Liuville problem (p8|), ( |3^ ) is determined functionally by the 
potentials in the bulk of the interval x G (—L/2, L/2). This makes it rather difficult to obtain the correlation functions 
by applying the well-elaborated methods of statistical analysis which are applicable to the evolution- type problems. 
One of the commonly used methods aimed at reducing boundary- value problems to evolutionary ones is the invariant 
imbedding method |2^, |27|. Yet in this study we apply a different approach which seems to be more general. We will 
search for Green function of equation ( p8| ) in the form 

Gii{x,x') =W^[i^+{x)^J-(x')e{x - x') +tp+{x')tp^{x)e{x' - x)] , (33) 

where ip±{x) are linearly independent solutions of the homogeneous equation (|28| ) supplemented with radiation con- 
ditions analogous to ( ^2|) at only one of the strip ends, x — ±L/2, in accordance with the "sign" index of ip±. The 
Wronskian of those functions is W, 9(x) is the Heaviside unit-step function. 

It is advantageous for the further analysis to represent the functions ip± (x) as superpositions of modulated harmonic 
waves propagating in opposite directions along the x-axis, 

V'±(2;) = 7r±(x) exp(±i/cix) — ij±{x) ex.p{^ikix) . (34) 

The WS approximation expressed in terms of inequalities ( |3l| ) suggests that the "amplitudes" 7r±(x) and 7±(x) in 
( ^ ) may be thought of as varying slowly as compared to the "fast" exponentials exp(±iA:ix). This makes it possible 
to obtain the "truncated" equations for those amplitudes by averaging the exact Schrodinger equation over "rapid" 



phases, as it is done in the theory of nonlinear oscillations [28 . Specifically, we multiply both sides of the homogeneous 



equation (^8|) from the left by the exponent function exp(=pifcix) and then average all the terms over space interval 
21, which can be chosen to have the arbitrary length between the micro- and macroscopic lengths of the problem, 

kp\r,^l^£sc,L . (35) 

The final result is not expected to depend exactly on the choice of the averaging interval. 

By means of such a sub-averaging, we arrive at the following set of first-order differential equations for the smooth 
amplitudes, 

±tt'^{x) + i'q{x)TT±{x) + C^{x)j±{x) ^ , 

(36) 

±7±(x) - iry(x)7±(x) + C±{x)tt±{x) = . 
The radiation conditions for ip± (x) are reformulated as the following "initial" conditions for the functions t:± and 7-1- , 

7r±(±L/2) = l, 7±(±L/2) = 0. (37) 
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STATISTICAL PROPERTIES OF SMOOTHED POTENTIALS 



At the last stage of developing the averaging procedure we specify statistical properties of the functions ri{x) and 
C±{x) entering equations (36). These random fields are defined as narrow packets of spatial harmonics of the initial 
potentials from the Ihs of Eq. (p8|). The contribution of intra-mode potential V{x) = Vh{x)+Vs{x) to these "smoothed" 
potentials is represented by the expressions 



rj^^^x) = 



1 

2h 



x+l 



x — l 



dt 



Vit) 



1 

2h 



x+l 



x — l 



dt 
2l' 



±2ikit 



Vit) 



(38) 



Unlike the local potential V{x), the operator potential JC, which is specified by the kernel (p9|), is, strictly speaking, 
non-local. Its part in the random fields r]{x) and C±i^) is described by more complicated expressions than those 
given by Eq. (|3^). The corresponding potentials, the same as the kernel ( [29| ) itself, are the differential operators. In 
calculation of different correlation functions the random fields ri{x) and C±i^) ^-re not of importance by themselves, 
but only their statistical moments. Subject to WS conditions, the latter can be calculated no matter what the locality 
of the corresponding potential may be. 

At weak scattering, all the potentials in Eq. ( |2^ ) may be thought of as Gaussian distributed Consequently, 
knowledge of binary correlators of those potentials is sufficient to govern the statistics of all physical quantities. It was 
shown in Ref. that under restrictions (^5|) for the averaging interval I, only a pair of correlators of the potentials 
r]{x) and C±{x) are not parametrically small, viz. {i]{x)r]{x')) and (^(±{x)(^^{x')). Regardless of the potential being 
local or not, calculation of these correlators yields 



{'q{x)r]{x')) = — F;(x - x') , 



(C±(:^)C;(^')) = ^^;(^-^') • 



The function 



(39a) 
(39b) 



Fiix) 



dq Jgx^^'^^ill) _ 1 \x\ 



2Tr [qlf 21 \ 21 



<d{2l-\x\) 



(40) 



in equations (^9|) is sharp in the scale of macroscopic lengths, so hereinafter it can be regarded as the ^-function in 
the "distributional" sense, Fi{x) d{x). 



The coefficients ^ in equations (39) represent the inverse scattering lengths, forward (/) and backward (6), 



respectively. They are contributed by all the potentials that result in both BH and BS scattering in equation ( |28| ) . In 
calculating these coefficients, it should be borne in mind that under WS conditions one can disregard the correlation 
betw een "height" potential (22a) and "slope" potentials (22b) and (p9[)[T^. This permits us to perceive the potential 
(22a), on the one hand, and potentials ( ^2b| ) and (29), on the other, as being associated with some different additive 
and mutually non-interfering scattering mechanisms. 

Scattering rates (both forward and backward) associated with "slope" potentials similar to (22b) and (^9|) were 
previously studied in Ref. [0 for CCB model of a roughly bounded strip. Although the CCB waveguide is somewhat 
different from the SSB waveguide addressed here, the difference appears to arise only in numerical factors at the 
potentials considered in paper and potentials (22b) and ( |29| ) of this work. Therefore, referring the reader to 
Ref. [|lOj for technical details, we present here the final expressions for the inverse scattering lengths associated with 
BS scattering, which are valid for the SSB waveguide model. 



-is) 



1 



2fc2 \D 



dq 
2^ 



~3 



2 J2 BL [(2fci + g)'aL")(fci +q) + (2fci - qfgl^Hh q) 



(41a) 



2k'f \D 



dq 
2^ 



[q' -kiYW{q-ki)W{q + ki) 



TT 

y 



m=2 



BUq' - kl)g^^\q) 



(41b) 
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The functions W{q) and ^^^(g) in Eqs. (gl|) are the Fourier transforms of the correlation function yV{x) from (^) and 

the evanescent Green function (^6|), respectively. 

As regards the BH scattering due to the potential (22a), the corresponding scattering rate is worth considering here 
in more detail. At first glance, it seems natural to find the corresponding frequency by expanding the potential ( 22a ), 
with condition (p^), in small fluctuations of the conductor width. 



Vhix) 



2^ 
D3 



Aw{x) 



(42) 



However, in trying to improve the obtained scattering frequency by retaining the terms of higher order in Aw{x) it 
turns out that the corresponding series converges non-uniformly, so that for its convergence it is necessary to hold 
a great number of terms. The similar problem was encountered earlier in Ref. ||30| , where the artificial "cutting" 
parameter was introduced for the corresponding series to become ultimately convergent. 

This difficulty can be overcome if the expansion of th e BH potential in series of small displacement of the confining 
surfaces is discarded, but the exact expression (22a) is used instead. This can be easily done for the gaussian 
roughness model, the appropriate scheme being presented in the Appendix. With this technique, one can obtain for 
inverse scattering lengths pertaining to the potential (22a) the expressions below, which are valid in the case of small 
boundary asperities, Eq. 



-{h) 



-(h) 



47r4a2 
47rV^ 



W(0) , 
W(2fci) 



(43a) 
(43b) 



It is noteworthy that extinction lengths (p3[), being obtained by rigorous calcul ation , coincide exactly in form with 
those obtained by means of the lowest-order expansion (^) of the Hamiltonian ( ^2a| ) in powers of /S.w{x). 



CONDUCTANCE AND RESISTIVITY MOMENTS 



To perform the statistical averaging based on the effective zero-scale correlation of random potentials (|39| ) it is 
necessary to express the conductance (^) in terms of smooth amplitudes 7r± and 7± . Under WS conditions (|3l|), when 
substituting the Green function in the form (^3|) into the expression ( |30| ) , it is sufhcient to differentiate over coordinate 
variables only fast exponentials in wave functions (34). As a result, the conductance takes on the intermediate form 



4kj 



L/2 



dx< 



-L/2 



[k+(x)p - |7+(x)p] / dx' [k-(x')P - 

J -L/2 

+ [k_(^)P - l7-(^)P] r'dx' [|7r+(x')P - \l+{x')\' 



(44) 



Within the same accuracy, the Wronskian W in (^4|) equals 

W « 2iki [7t+{x)tt_{x) + -f+{x)'y-{x)] 2ifci7r±(=Fi/2) , (45) 

where the last equality is a consequence of "boundary" conditions ( |37| ) , which are valid within WS approximation as 
well. 

Formula ( ^ ) can be simplified if the symmetry properties of equation ( ^ ) and hermicity of the corresponding 
ID Hamiltonian are taken into account. Since the problem (|3^), ( |37| ) is of evolutionary type, its solution can be 
represented in terms of an a;-ordered matrix exponential. 



I± (x) = exp 



± 



±L/2 



dx'h{x') 



Here the matrices of smooth amplitudes, I±(a;), have the form 



I+(x) 



7+ (a;) 7r;(a;) 



I_(.x) = 



7r_(x) 7l(a;) 
7_(a;) 7rl(x) 



(46) 



(47) 
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h(x) is the random field matrix, 



b(x) 



ir]{x) C+{x) 
C-(x) -ii]{x) 



(48) 



whose ofF-diagonal elements are interconnected by the equality C-{x) = C+{^)- The operator Tx in Eq. ([iq ) arranges 
the multipliers in each of the terms of the exponential series in order of decreasing their coordinate arguments from 
left to right. 

Taking advantage of the operator identity Indet A = Spin A, matrices (E^) can be shown to be unimodular, 



This relation, along with equality 



detl±(x) = |7r±(x)p - h±{x)\^ = 1 . (49) 
|), results in the following form of the conductance of a single- mode wire, 



giL) = \7rg\TL/2)\ 



(50) 



Holding to Landauer's concept, the quantity n^^{^L/2) is to be interpreted as the transmission coefficient of a single- 
mode quantum waveguide of length L. This interpretation is supported by the following argumentation. From the 
structure of wave functions (^) it follows that the ratio r+ (x) = 7+ (a;) /7r_|_ (x) is defined as the reflection coefficient for 
the harmonics ki incident onto the interval (x, L/2) with a unit amplitude from the left-hand side. Correspondingly, 
T-{x) = 7_ (a;)/7r_ (a;) represents the reflection coefficient of the harmonics —ki incident onto the interval {—L/2,x) 
from the right-hand side. Subject to this definition, equation (^9|) can be re-written in the form of the conservation 
law in a non-dissipative medium, 



|r±(x) 



1 



(51) 



whereupon the interpretation of quantity (|50|) as a square modulus of the transmission coefficient of the disordered 
interval {—L/2, L/2) seems to be apparent. 

With the conservation law (^l|), it is convenient to perform subsequent calculations of the statistical moments of the 
conductance using the equation for r-|-(a;) rather than that for the transmission coefficient 7r^^(a;). From Eqs. (|3^), 
the reflection coefficient can be found to obey the Riccati-type closed equation subject to the zero initial condition. 



± 



dT±{x) 



dx 

r±(±L/2) 



= 2i7j{x)r±ix) + C±ix)Tlix) - C±(a;) 
. 



(52) 



The forward-scattering random field r]{x) may be eliminated from Eq. ( p^ ) by concurrent phase transformation of the 
function r±(x) and the backscattering field C±{x), 



r±(x) = r^'™>(x)exp 



C±(a;) = C±"''"'^(a;)exp 



±2i 



±2i 



±L/2 



dx'ri{x') 



dx'ri{x') 



±L/2 



(53a) 
(53b) 



This transformation keeps the conductance (|5C|) and correlation relation ( 39b ) unaffected, so one may put function 
77(2;) in Eq. ( [s^ ) equal to zero. As a consequence, the outcome for an arbitrary moment of the conductance is to be 
specified exclusively by the backscattering of the electrons, i.e. it will depend on the scattering length which inverse 
value is the sum of inverse lengths ( [41b| ) and (43b). 

To proceed further, consider the nth moment of the local reflection coefficient squared modulus, 

R^{x) = {\r±{x)n . (54) 

Since the stochastic problem ( p2[ ) is of the evolutionary type, it can be reduced, via Furutsu-Novikov formalism [p7[, 
to the differential-difference equation for moments (p^) (see also Ref. [pl|). 



± 



dx 



[Rt+Ax)-2Rt{x)+R^_,{a 



0,1,2, 



(55) 
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Here L^^^ = L^f' ^ + l'^^^^ ^. The initial condition on the coordinate x to the solution of Eq. ( p5| ) is 

i?±(±L/2) = (5„o. (56) 



As for the dependence of R^ix) on the discrete variable n, it follows from the definition (54) that Rq{x) = 1 and 
Rn{x) ^ as n ^ oo. 

The solution of Eq. (js^) that matches all the above mentioned conditions can be expressed through the probability 
function P^{u\x) and represented, upon due parametrization, in the form 

Rtix)^ / duPi{u\x)(-—) . (57) 



,"" + 1, 

Correspondingly, the statistical moments of the conductance ( |50| ) are represented by the following integral, 

(5"(i)) = ( (1 - |r±(TL/2)H" ) = duP^{u\ T m (^)" ■ (58) 

To obtain the probability density P^{u\x) one should substitute R^{x) in the form ( |57| ) into equation (|55|), thus 
obtaining the Fokker-Plank equation, 

dP^ = _ _ ,.dP^ 

" dx du^ ' du ' ^ > 

which is supplemented, according to Eq. (|5^), by the initial condition on the coordinate a;, 

P^{u\±L/2) = 5{u-l-Q) . (60) 

Besides, normalization of the function P^{u\x) to unity is ensured by R^{x) = 1. This implies that the distribution 
function is integrable over the variable u, in particular, at u ^ 1 and u ^ oo. 

The solution to Eq. (p9|), which meets the above mentioned requirements, can be found by using the Mehler-Fock 
transformation and has the conventional form IBSI 



,± , , , 1 fL^2x\ f Lt2x\ [°° vdv 



P,(cosha|.) = ^(^^j expj^-^^jy^ (^^ih^^r^^ih^ 



4 V 2Lh 



(61) 



where the change of a variable has been made u — cosh a, a > 0. With this expression, equality ( |5^ ) yields the 
relatively simple (as well as suitable to analyze) formula for the nth moment of the dimensionless conductance. 



= ^(fj "H-3i;jy. ij/.*'™'' '''' 

n = 0,±1,±2,... . 

The result ( |62| ) completely determines main averaged transport characteristics of a single-mode conducting strip. 
In particular, although the conductance itself is not a self-averaged quantity, one can, in principle, calculate its 
self-averaged logarithm using the whole set of statistical moments (^2|). However, it is much easier to obtaine this 
quantity directly from equations (^6|), omitting cumbersome manipulations with a logarithmic series of the terms (|6^). 
Specifically, by differentiating the quantity |7r^^(a;)p over x we arrive at the equation 

±^ln|7r±i(x)p=Ci(x)r±(x)+C±(x)ri(x) . (63) 



By integrating Eq. ( |63| ) over the interval (— i/2,L/2), with ( |37| ) taken into account, wc immediately obtain the 
logarithm of the conductance in the left-hand side. Before averaging the terms in the rhs of Eq. (|63|), we point out 
that from equation ( |52| ) it follows that function T± (x) , being considered as a functional of random fields C± and C± , 
depends on the value of these fields exactly within the interval (x,±L/2), according to the sign index (±). Since at 
weak scattering all the above mentioned random fields can be regarded as Gaussian distributed functional variables. 
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we can apply the Furutsu-Novikov formalism for its averaging! 27 . The average of the first term in the rhs of Eq. (|6 
can, therefore, be presented in the form 



(Ci(x)r±(x)) = ± [^'^'\x' {cUx)c±{x')) 

J X 



dT±{x) \ ^ J_ / dV±{x) 
dC±{x')/ 2Lb\dC±ix') 



(64) 



where we have used the effective (5-corrclation of the fields C± (x) and (x) . The variational derivative in (^J) can be 
readily calculated with the use of equation (^2|), and thus it turns out to be equal exactly to unity. Finally, we arrive 
at the well-known result for ID disordered systems, 



{lng{L)) = -L/Lb , 



(65) 



signaling the exponential fall of the conductance with a growing length L at the so-called representative (non-resonant) 
realizations of the random potential |29]. 



DISCUSSION OF THE RESULTS 



With general formula (^2|), write down the expressions for average resistance {g ^{L)) and average conductance 
{g{L)). At n = — 1 the integrals in Eq. ([6^) can be calculated exactly, so the average resistance is equal to 



{9-\L)) = i 



1 



2L 



(66) 



At n = 1 the integration can be performed asymptotically in the parameter L/Lb, giving rise to the following 
expression: 



{9{L)) 



1 - L/Lf, if L/Lf, < 1 

2-i7r5/2 [L/LbY'^'^ exp {~L/ALb) if L/U > 1 



(67) 



The results (p5D~(|67[) are completely in line with the concepts of the localization theory for one-dimensional dis- 
ordered systems. Obvious indications of the ballistic electron transport in short wires can be easily seen, g{L) « 1 
at L <C Lb- Also, no signs of diffusive motion of the electrons in long wires are present in any result. And con- 
versely, in long wires, L ^ Lb, the resistance ( |66[ ) displays an exponential increase with a growing strip length, and 
the asymptotic ( |67| ) shows an exponential decrease of the average conductance as the length L exceeds the value of 
ALb- That behaviour is characteristic for conduction electrons undergoing Anderson localization. The inverse of the 
quadruple backscattering extinction length is equal to the Lyapunov exponent for the electron wave function in the 
case of dimension unity ||29| , so that the quantity £\oc = 4Lf, is conventionally called the (one-dimensional) localization 
length. 

Both of the formulas (|6^) and (|6^), and, hence, (|6^ ) and (|67|), are universal in that they are applicable for any 
one-dimensional degenerate system subject to weak static disorder. Only the particular dependence of scattering 
lengths Lf^b and, consequently, the localization length £ioc is determined by the physical nature of disorder. In this 
paper, we have found that if the boundary roughness prove to be just the main cause of the disorder in a 2D single- 
mode conducting strip, the interpretation of scattering mechanism may be substantially different, depending on the 
interrelation between by-slope scattering lengths (|4^) and by-height lengths (43). To make the correct comparison 
of the lengths the roughness statistics needs to be specified. We will make the comparison for two characteristic 
models consistent with the demand for analyticity of the random functions ^1,2(2^) descriptive of boundaries of the 
conductor. Specifically, we examine the case of asperities subject to gaussian (exponential) correlation statistics. 



W(a;) = exp (— a;^/2r^), an those described by the Lorentz (power-type) correlation function W{x) = [l + [x/tcY 



For the gaussian roughness, from (41b) and 



-{h) 



Gauss 



D2 




the estimates can be obtained 



-2fc^2A 



if Tc/D < 1 (i.e. kirc < 1) 
if rc/D > 1 (i.e. kirc > 1) 



(68a) 



■is) 



Gauss 



4 



1 (fciTc)'' exp (— fcj? 



if rc/D < 1 (i.e. kirc < 1) 
if r^D > 1 (i.e. kirc > 1) 



(68b) 
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If the hight correlation is lorentzian, the estimates change to the following, 



-(h) 



Lorcntz 



1 



D) 



exp 2fcirc 



if Tc/D < 1 (i.e. kirc < 1) 
if Tc/D > 1 (i.e. kirc > 1) , 



(69a) 



Lorcntz 



4 

(7 \ Tr 



(fcirc)^exp (-2fcirc) 



if rc/D <C 1 (i.e. kirc < 1) 
if rc/D > 1 (i.e. kirc > 1) 



(69b) 



Based on estimates (|68|) and (|69|), the relative intensity of BS and BH scattering can be estimated as follows, 



L 



(h) 



Gauss 



f- 



d) \r,. 



if rc/D < 1 (i.e. kirc < 1) 
if rc/D > 1 (i.e. kirc > 1) , 



(70a) 



-(h) 



Li 



Lorentz 



if rc/D < 1 (i.e. kirc < 1) 
if rc/D > 1 (i.e. kirc > 1) 



(70b) 



From these estimations it can be seen that in all of the limiting cases considered here the characteristic ratio L^'*' /i^"' 
is determined by the product of small Rayleigh parameter [a/D^ and some scale parameter, though individual for 
different roughness correlation, which depends on the relation between the roughness correlation length and de-Broglie 
wave length of the electrons. In the case of small-scale asperities, when kprc ^ 1 (or, which is the same, rc/D <C 1), 
regardless of the correlation model the relative intensity of BH and BS scattering is characterized by the parameter 
(cr/_D)^ [D/rc)'^ = (cr/rc)* / {(j/ D)^ which varies over a rather wide range. This is because of BH and BS scattering 
being associated with independent physical origins. Whereas scattering from the potential Vh{x) is governed by the 
height of the boundary roughness and thus is estimated mostly in terms of the parameter cr/D, Eqs. (38a) and (69a), 
the BS scattering is mostly determined by the slope of the asperities, i.e. by gradients 2(2^)1 ^^^id is, consequently, 
governed by the parameter a/rc, Eqs. (68b) and (69b). 

It would be tiresome to discuss here in detail the interrelation between BH and BS scattering mechanisms assuming 
the asperities to be large-scale, when the parameter kprc ^ 1 (i.e. rc/D 3> 1), since in this case the the result 
depends largely on the correlation model. We leave this particular analysis to an interested reader. 

It is noteworthy that in the event when the "slope" mechanism dominates the "height" mechanism, the "surface" 



scattering rate is 
diffraction theory 



Droportional to the forth power of the rms height a rather than to cr^, which is customary in the 
T^ . This fact must be taken into account when analyzing experiments aimed at reproducing the 
surface shape using the data on quantum, as well as classical, wave scattering in rough-bounded waveguide systems. 

Apart from the identification of dominant scattering mechanism, specifying scattering lengths ( |68| ) and (|69| ) for 
different statistical models of boundary roughness also allows for the criteria of validity of the obtained results in 
terms of essential physical parameters of the disordered systems. Along with the presumption of smallness and 
smoothness of boundary asperities (see Eqs. ( |l^ ) and (^9|)) the criteria are dictated by ( ^l|) of weak scattering. In 
specifying those criteria, the smallest of the extinction lengths L^*'*^ should be taken to substitute the length £sc m 
( ^ since the inequality always holds true Lf < Li,. 

In conclusion, we make some remarks concerning the methodological side of the problem of wave scattering from 
rough waveguide surfaces. As far as we know, until the present time there has not been made any reasonable distinction 
between BH and BS scattering in such systems. Only the existence in general of different competing mechanisms 
responsible for wave scattering from rough surfaces was indicated in Ref. ]T6t on the basis of the experimental 
results. Accordingly, the relative function of these scattering mechanisms in dynamic processes in waveguide-like 
systems was not properly analyzed. Meanwhile, in the course of this work we have made certain that application of 
linearized (impedance-type) boundary conditions to a single-mode waveguide is equivalent to retaining in Eq. ( p^ ) 
the approximate potential (^2|) instead of its exact value ( ^2aD , and also disregarding all "slope" potentials. However, 
omitting the latter potentials implies the neglect of BS scattering mechanism, that is proven to be not always justifiable. 
Alternative small-slope approximation of Refs. ||l^, being guided by a solely slope parameter of the rough surface, 
does not allow one to separate BH and BS scattering mechanisms as well. 

In this work, the method has been suggested within the framework of which both of the above-mentioned scattering 
mechanisms appear quite naturally, being associated with different terms of the Hamiltonian. We have demonstrated 
that at least in a single-mode waveguide the scattering caused even by mildly sloping boundary asperities can be 
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attributed to either "height" or "slope" scattering mechanism, depending on the statistical properties of the roughness. 
The competition between these mechanisms is governed by physically different parameters. It is noteworthy that taking 
into account the BS mechanism is particularly essential if boundary asperities are classified as being large-scale. 
N.M.M. acknowledges support from CONACYT. 



APPENDIX: CORRELATION CHARACTERISTICS OF BY-HEIGHT SCATTERING 



The effective sub-averaged potentials rjh{x) and C,h±{x) corresponding to the original potential ( P2a|) have the form 



TT^ dt 



1 



1 



To calculate the averages in Eqs. (71), it is convenient to make use of the method below, 

1 



d 



1 



w{t) 



1 \ 



1 



dD \ w{t) 
i dq (^exp {iq[D + A£,{t) + iO]} 



(71a) 
(71b) 

(72) 
(73) 



In the case of boundary roughness obeying gaussian statistics, the integral in (|7^) can be averaged without difficulty. 
Since for the SSB waveguide the equality holds true 



the correlator (73) equals 



At) 



2a y 2 



exp 



_D2 



erfc 



-iD 



(74) 



(75) 



where erfc(...) is the probability integral (see, e.g., Ref. p2[). 

With the use of ([75|), the correlator (|7|) can be accurately computed. However, subsequent analysis can be 
performed analytically only in the case of small- height roughness. With inequality (p!^), the averages (^) and ( [T^ ) 
are asymptotically equal to 



1 



wit) 
1 

w^t) 



1 

D 

1 



1 + 4 



D 



(76) 
(77) 



The next thing to proceed is the calculation of binary correlation functions {rih{x\)'r\h{xi)) and {C^h±{x\)Qi±k^'2S) ■ 
In doing so, one has to calculate the correlator 



£(ii - ^2 



1 



1 



W2(ii) \vflit{) 
1 1 



U.2(t2) \zi;2(t2) 



l«2(ti)u,2(t2 



W2(fi 



W2(t2) 



(78) 



In the case of small-amplitude roughness, the second term in the rhs of Eq. ( |78| ) is governed by asymptotic (77). As 
for the first term, for its calculation the method which was already applied to calculate the average (^2|) is helpful. 
To this end, represent the desired correlator in the form 



1 



1 



,«;2(ii)M;2(t2 



a2 



lim „ „ „ 
D'~.D dD dD 



- jj^ dqidq2 (exp [iqi [D + A^{ti) + lO] - iq2 [D' + AC{t2) ^ iO]}) 



(79) 
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Averaging of the exponent function in (|79|) yields readily 

(exp [^9lA^(^l) - iq2A^{t2)]) = exp {-2 {qj + ql) + 4<?ig2(T2>V(ti - <2)} , 
whereupon we arrive, after some tedious manipulations, at the integral representation 



1 



1 



1 



{t\)w'^{hl {2gY^W+W- Jo 



du cos 



D 



dv 



ip{u) 



V U 



exp 



Here, for the sake of the formula compactification, we use the notations 

(p{u) = u^/W-/W+ , W±^l± Wih - <2) 



(80) 



(81) 



(82) 



Correlation function (|8l| ) can be calculated numerically for the arbitrary value of a/D. But in the case of small- 
height roughness the integration in (81) can be performed analytically, yielding the asymptotic result 



1 



1 



By substituting (||) and (|7^) into (^ we obtain 



1 



16cr^ 



W(tl-t2) 



(83) 



(84) 



After simple integration of correlator (^J), resulting from definition (|71|), we arrive at the final expressions for the 
required correlators. 



1 



{Vh{xi)llh{x2)) = -^Fl{xi - X2) 



{C,h±{xi)Ch±{x2)) 



-(h) 



Fi{xi - X2) 



(85a) 
(85b) 



In Eqs. (|8^), the extinction lengths L^l^\ associated with the height potential (22a) are given by Eqs. ([43|). 
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